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Abstract 

We give a formulation of quantum ergodicity for Pauli Hamiltonians with arbi- 
trary spin in terms of a Wigner-Weyl calculus. The corresponding classical phase 
space is the direct product of the phase space of the translational degrees of freedom 
and the two-sphere. On this product space we introduce a combination of the trans- 
lational motion and classical spin precession. We prove quantum ergodicity under 
the condition that this product flow is ergodic. 
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1 Introduction 



Quantum ergodicity describes the equidistribution of eigenfunctions in the semiclassical 
limit (h — > 0), in the sense that almost all expectation values of quantum observables 
tend to a classical mean value of the corresponding classical observable. In particular, 
the Wigner functions associated with almost all eigenfunctions of the Hamiltonian become 
equidistributed on the energy shell in phase space. A sufficient condition for quantum 
ergodicity to hold is ergodicity of the corresponding classical dynamics, as was first stated 
by Shnirelman [|TJ. It is remarkable that quantum ergodicity is one of the few results in 
quantum chaos for which there exist mathematical proofs; the first complete ones were 
given by Zelditch |J and Colin de Verdiere ||. These are results for scalar Hamiltonians, 
i.e. for the Laplacian on compact manifolds or, more generally, for Schrodinger operators 
in IR d ||]. The case of quantum billiards was considered in [EL EL 

If one wants to describe particles with internal degrees of freedom, such as spin, one 
has to deal with matrix valued Hamiltonians. Since spin is a purely quantum mechanical 
property, it is a priori not obvious what should serve as the corresponding classical sys- 
tem. In |7j] quantum ergodicity for Pauli Hamiltonians with spin 1/2 was proven under 
the condition that a combination of the classical translational dynamics and a quantum 
spin dynamics driven by the translational motion is ergodic. This kind of mixed clas- 
sical/quantum description in terms of a skew product dynamics arises naturally in the 
semiclassical analysis of Hamiltonians with spin || |9], |10| . 

Although spin is a genuinely quantum mechanical property, there exists an intuitive 
classical analogue given by the so-called vector model which provides a picture of spin as 
a vector of fixed length whose dynamics is similar to that of angular momentum. The 



mathematical description of this model goes back to Thomas [11 1, who gave an equation 



of motion for a classical particle with a magnetic moment precessing in external fields. In 
a semiclassical analysis of spinning particles various objects can be related to this classical 
spin precession [[12|, ||, [L3[]. In the vector model the equations of motion are of first order in 
time, and for pure spin precession can be formulated in a Hamiltonian framework. Thus, 
after suitable normalisation of the spin vector, the unit sphere S 2 serves as the phase space 
of classical spin precession. 

Based on ideas of Stratonovich fH |, Gracia-Bondfa and Varilly []15 , |16f developed a 
Wigner- Weyl calculus for arbitrary spin j, in which quantum mechanical observables, i.e. 
hermitian (2j + 1) x (2j + 1) matrices, are represented in terms of functions on the phase 
space S 2 . We use this description to rephrase the problem of quantum ergodicity for non- 
relativistic particles with spin and to prove for arbitrary spin the respective statement 
under the condition that a different skew product of translational and spin dynamics, 
where now the spin part is given by classical spin precession, is ergodic. 

The paper is organised as follows. In section |2] we briefly review properties of the Weyl 
correspondence for systems with only either translational or spin degrees of freedom. In 
section ^ we then discuss how in the semiclassical limit of systems with both kinds of 
degrees of freedom the motion is governed by the skew product dynamics. We also state a 
corresponding Egorov theorem. Section |] is devoted to the proof of quantum ergodicity for 
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Pauli Hamiltonians with arbitrary spin. Some particular representations of the Wigner- 
Weyl transform for spinors and a relation between ergodic properties of the two types of 
skew products are discussed in two appendices. 



2 Weyl quantisation and classical limit 

The Weyl quantisation of systems in M. d without spin, i.e. taking only translational de- 
grees of freedom into account, is based on unitary irreducible representations of the Weyl- 
Heisenberg operators 

p(p, q) := eK*****) , (2.1) 

where Pk, X^, k = 1, . . . , d, are momentum and position operators, respectively, satisfying 
canonical commutation relations. According to the Stone-von Neumann theorem such 
representations are unitarily equivalent to the Schrodinger representation ps(p, q) of (|2.1| ) 
on the Hilbert space 7itrans — L 2 (M. d ), in which X/~ is realised as a multiplication operator 
and Pk = faf - - The Schrodinger representation hence acts on ip e L 2 (IR d ) as 

(ps(p, q) ij>) (*) = efr^to ^{x + q) . (2.2) 

See |L7j for further details. 

A classical observable B(p, x) is a real valued function on phase space M. d x M d , to 
which the Weyl quantisation assigns a symmetric operator B on L 2 (M. d ), 

B-T^Tuf I B(p,q)p s (p,q)dpdq, (2.3) 

where 

B(p, q) ■= tt^-j / I B{£, z) e-**^*") d£ dz (2.4) 
{2nh) d J Rd J Rd 

denotes the Fourier transform of the classical observable. The operation of B on (suitable) 
•0 G L 2 (R d ) then reads 

B^(x) = -±— d f [ e^(*-y) B (p,i(x + y))iP(y)dydp, (2.5) 
(Ann) j Rd j Rd 

In fact, B can even be a (tempered) distribution on phase space, B e S'(M. d x M d ); in this 
case the operator B is defined on the domain 5(R d ) C L 2 (M. d ). Reversing the above rea- 
soning, one notices that the Schwartz kernel theorem allows to represent every continuous 
linear map B from S(R d ) into S'(R d ) Weyl operator (|2.5|) with symbol B(p, x), which 
in general is a distribution on phase space (see, e.g., |T7| ). The map B \— > B(p,x) can be 



3 



made more explicit if B is a trace class operator on L 2 {R d ). Since ps(p,o) is bounded, 
B PsiPi Q) is also of trace class, and 

B( P ,q)=Tr(Bp s (p,q)) , (2.6) 

where Tr(-) denotes the operator trace on H trans . We remark that the symbol of a Weyl 
operator can in principle depend on h. In such direct interpretation of B(p, x) 

as a classical observable is inappropriate. However, if one restricts the class of observables 
suitably, the symbols possess asymptotic expansions as h — > 0, 

B(p,x;h)~J2ti k Bk(p,x) ; (2.7) 

fc>0 



for the precise meaning of (|2.7|) see, e.g., Now the principal symbol B (p,x) 

is a function on phase space independent of h that serves as the corresponding classical 
observable. 

Summarising its most important properties, the Weyl quantisation can be characterised 
by the following properties (see, e.g., [[H], |T^, [TP|]): 



(i) The map B i— > B(p, x) is linear. 



(ii) £t ^ B(p,x). 

(iii) If i? is of trace class, then 

TtB = — i— / / B(p,aj) dp da; . (2.8) 
(27r^) d y R<i J R d 

(iv) Covariance: If G G Sp(2cf, R) is a linear canonical transformation with metaplectic 
representation n M (G), then B \— > B(p,x) implies 

n M (G^Bn M {Gf) ^ B(G(p,x)) . (2.9) 

We now turn attention to a quantum mechanical spin, without translational degrees 
of freedom. Its (pure) states are described by vectors in the Hilbert space TC sp in = C 2j+1 , 
where j G |N denotes the (fixed) spin quantum number. Spin observables are therefore 
hermitian (2j + 1) x (2j + 1) matrices. The corresponding classical phase space is given 
by the two-sphere S 2 , considered as a symplectic manifold whose symplectic structure is 
provided by the area two-form. In spherical coordinates (9, <p) for S 2 a symplectic chart is 
given by (p,q) = (—cos9,(f) G (—1,1) x (0,27r). For reasons of convenient normalisation 
we here choose the symplectic two-form as 

^G°>9) := r =dpAdg= 1 sin^d^Ad^ . (2.10) 
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In the following we will often denote points of S 2 by n G R 3 with \n\ = 1 so that in spherical 
coordinates n = (sin#cos<^, sin # sin y?, cos#) T . In these coordinates the normalised area 
two-form is given by d/i S 2(n) = ^sin^d^ A dtp. Below we will show that the quantum 
mechanical spin operator S is then associated with the vector s := hWj(j + 1) n G M 3 of 
length hyj j(j + 1); in other words, n = s/\s\ G S 2 represents a normalised 'classical' spin 
vector. 

A Weyl correspondence for spin, with properties as close as possible to the ones dis- 
cussed above in the case of translational degrees of freedom, has to assign a real valued 
function on S 2 to a hermitian (2j + 1) x (2j + 1) matrix in such a way that properties 
analogous to (i)-(iv) of above hold. Based on ideas of Stratonovich [14]], Varilly and Gracia- 



Bondia |15], |TB| have constructed such a Weyl correspondence. In their approach quantum 
mechanical observables b are represented as 

b = (2j + 1) / 6(n) A» d// s2 (n) . (2.11) 
7 s 2 

Here Aj(n) is a function on S 2 taking values in the hermitian (2j + 1) x (2j + 1) matrices, 
and the symbol b(n) G L 2 (S 2 ) can be obtained from the observable b in analogy to 



6(n) =tr(SA i (n)) , (2.12) 

with tr(-) meaning the matrix trace. In jl6| it is shown that there are 2 2j possibilities for 
a Weyl correspondence with a kernel Aj-(n) that fulfills 

(a) Aj(n)1 = Aj(n) for all n G S 2 . 

(b) (2j + 1) / ga Aj(n) d / u s2 (n) = l 2i+1 . 

(c) Covariance: If g G SU(2) and 7Tj denotes the (2j + l)-dimensional unitary irreducible 
representation of SU(2), then Aj(ip(g)n) = ^j(g) Aj(n) iTj(g)^- 

In the covariance property of the kernel (f denotes the covering map from SU(2) to SO (3): 
For every g G SU(2) one defines the adjoint map Ad 9 as the linear map of the Lie algebra 
su(2) given by X i— > g' X g. Expanding X G su(2) in terms of the Pauli matrices, X = x cr, 
x G M 3 , the adjoint map operates as Ad g (x ■ er) = (ip(g)x) ■ cr. In this way every g G SU(2) 
is mapped to some (p(g) G SO (3); in fact this map is two-to-one and provides the universal 
covering of SO(3) by SU(2). Hence (p(g)n means the rotation of n G M 3 (with |n| = 1) by 
ip(g) G SO(3). In |16| one choice out of the 2 2j possibilities to define Aj(n) is singled out 
due to its connection with spin-coherent states; however, for our purposes this particular 
choice is not important and, thus, will not be specified here. For illustration, some explicit 
examples are shown in appendix [A] 

The properties (a)-(c) of Aj(n) listed above now imply properties of the symbols b(n) 
that are closely analogous to those of symbols in the case of translational degrees of freedom: 
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(i) The map b i— > bin) is linear. 

(ii) i — > b(n). 

(iii) trS = (2j + 1) / s2 6(n) d/i S2 (n). 

(iv) For every g G SU(2): 7r i (g') t bir^g) i-> b(ip(g)n). 
In addition one obtains the simple relation 

tr(aS) = (2j + 1) / o(n)6(n) d// S a(n) (2.13) 
is 2 

for the trace of a matrix product. 

The covariance of the quantisation as expressed by (iv) relates the quantum mechanical 
time evolution of a spin to a classical dynamics (see also [[RJ): Let diTj denote the (2j + 1)- 
dimensional derived representation of su(2), i.e. diij(X) = lj^iTj(e lXX )\\ =0 . Then the 
spin operators Sk, k = 1,2,3, are given by Sk = |d7Tj((T fc ), with commutation relations 
[Sk, Si] = ihSkimSm- The Weyl correspondence now assigns the vector valued symbol s = 
hyJQ + 1) n to the vector S of these quantum mechanical spin operators, see appendix 

E 

A typical (time independent) quantum Hamiltonian is given by 

H spin = S-C = ^C-dn J (a) , (2.14) 

where C G K 3 is some constant vector. This Hamiltonian generates a quantum mechanical 
time evolution that is governed by the unitary operators 

e -iV' = ff3 (e^ c,Tt ), (2.15) 

so that a time evolved observable reads 

b{t) = e ^ int be-^ int = ir j (e-i c -' Tt ybir j (e-i c -' Tt ) . (2.16) 

If b(n) is the classical observable associated with b, the relation (|2.16|) and the covariance 
property (iv) hence assign the classical observable 

b{t){n) = b((p(e-* c -' rt )n) (2.17) 

to b(t). In this way one is provided with a dynamics n i— > n(t) = {p(e~^ c ' eTt )n on S 2 that 
is governed by the equations of motion 

h(t) = Cxn(t) (2.18) 

describing a precession of n{t) about C. On the other hand, the Weyl correspondence 
associates the classical Hamiltonian H sp - m (n) = s ■ C to the quantum Hamiltonian ( |2.14| ). 
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Together with the Poisson bracket {•, -}s 2 derived from the symplectic structure on S 2 given 
by uj one therefore identifies (|2.18|) as the equations of motion generated by H sp i ni i.e. 



h{t) = {H spin , n(t)} s2 = C x n(t) . (2.19) 



Thus the assignment (|2.12|) of a function on S 2 to a quantum observable commutes with 



the (quantum mechanical or classical, respectively) time evolution, see also appendix [B]. 
This property is a special feature of the Weyl correspondence for spin. In the case of 
translational degrees of freedom the analogous relation, known as the Egorov theorem (^Of , 
only holds for the assignment of the principal symbol B (p, x), and not of the full symbol 
B(p,x), to a quantum observable. 



3 Coupling of spin and translational motion 

The quantum mechanical description of particles with spin requires to couple translational 
and spin degrees of freedom. The relevant Hilbert space 7i therefore has to be the tensor 
product of 7Y t rans and 7"4 P in, i-e. Ti = L 2 (IR <i )®C 2j,+1 . The state vectors of a spinning particle 
are hence (2j + l)-component spinors whose components are square integrable over M. d . 
Observables are represented by self-adjoint operators on Ti, which in Weyl representation 
read 

^ x ) = 7^dl I e^<^B(p,l(x + y))^(y)dydp. (3.1) 
[ZTin) j R d j R d 

Here the symbol B(p, x) is in general a (2j + 1) x (2j + 1) matrix valued distribution on 
phase space R d x R d . The operators (|0|) then act on ip e S(R d ) C 2j+1 . 

In a non-relativistic context the quantum dynamics is generated by the (Pauli) Hamil- 
tonian 

h 

Hp = H trans l 2j+1 + - dTT, (<x) ■ C (3.2) 

that consists of a scalar part H trans describing the dynamics of the translational degrees 
of freedom, and a genuinely matrix valued part that contains a coupling of spin to the 
translational motion; here Ck are Weyl quantisations of suitable functions Ck{p,x) on 
phase space. Typically, H trans is of the form 

4rans = ^~ (~ V - -A{xj) ' + e<f)(x) , (3.3) 

2m V 1 c / 

with (static) electromagnetic potentials A, <p, i.e. it is a Weyl quantisation of the classical 
Hamiltonian 
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Examples for the coupling of spin to the translational motion are given by a coupling 
to an external magnetic field, C(p,x) = —-^—B(x), or a spin-orbit coupling C(p,x) = 
4m 2 c 2 |a;| & ti\x\^ ( x x P)- ^ or ^ ne following the specific form of the Hamiltonian, however, is 
not relevant. We only require H (p,x) to fulfill certain criteria that guarantee Hp to be 
(essentially) self-adjoint, bounded from below, and to possess a purely discrete spectrum 
in some interval [E — s, E + s] . For details see 0]. 

For simplicity the further observables that we are going to consider shall be bounded 
Weyl operators on 7i. Their symbols are smooth, matrix valued functions on phase space, 
which may also depend on H in such a way that an asymptotic expansion (h — > 0) 

B(p,x;h)~ y %2?i k B k (p,x) (3.5) 

fc>0 

analogous to ( |2.7| ) holds (see for further details). As far as the translational degrees 
of freedom are concerned, the classical observable corresponding to B is then given by its 
principal symbol B (p, x), which is a function on phase space taking values in the hermitian 
(2j + 1) x (2j + 1) matrices. Thus spin is still described on a quantum mechanical level. 
As an example of an observable, although an unbounded one, that illustrates this situation 
consider the total angular momentum J = L + S, whose symbol 

J{p,x) = (x x p) t 2j+1 + h-dn j (a) (3.6) 

has a principal part that consists of orbital angular momentum, and a sub-principal part 
given by spin. 

The relation between the quantum and classical time evolution, i.e. the relevant Egorov 
theorem, has in this context been derived in [J7[|. In order to state this version of the Egorov 
theorem, let d(p, x, t) e SU(2) be the solution of the spin-transport equation 

d(p,x,t) + ^C($ t (p,x))-crd(p,x,t) = with d(p, x, 0) = 1 2 , (3.7) 

where x) = (p(t),x(t)) denotes the classical flow on the phase space M. d x M d gen- 
erated by the classical translational Hamiltonian H (p,x), i.e. (p(t),x(t)) are solutions 
of Hamilton's equations of motion with initial condition (p(0),x(0)) = (p,x). Hence, a 



comparison with (|2.15| ) reveals that 7Tj(d(p,x 1 t)) is the quantum mechanical propagator 
for the spin degrees of freedom along the trajectories $*(p, x) of the classical translational 
motion. The Egorov theorem relevant in the present situation now states that the time 
evolution B(t), generated by a Pauli Hamiltonian ( |3.2|) , of an observable ( |3.1|) with a sym- 
bol allowing for an asymptotic expansion (|3.5|), is again an operator of the same type and 



its principal symbol B (t)(p } x) reads 

B (t)(p,x) = -K,{d{p,x,t)) ] B^ip.x^-K^dip.x.t)) . (3.8) 

In this expression the translational degrees of freedom are obviously propagated classically 
by means of the flow $*. In contrast, a glance at the relation ( J2.16 ) reveals that the 



8 



dynamics of the spin degrees of freedom are genuinely quantum mechanical, but are driven 
by the classical translational flow and take place along the trajectories of this flow. This 
observation reflects the mixed (classical/quantum) level of description of the translational 
and spin degrees of freedom. 

We now proceed to describe spin within the framework of the Weyl correspondence 
outlined in the previous section. Thus, according to ( |3.1| ) and ( 2.1 1|) quantum mechanical 
observables B are represented as 

^ x ) = TTmf I I ^ pix - y H(p, 1 1 (x + y),n)A 3 (n)^(y)dydpdfi s2 (n), (3.9) 
(2nh) a J Rd J R d J S 2 

where according to (|2.12| ) 

b(p, x, n) = tv(B(p, x) Aj(n)) (3.10) 

is the associated scalar symbol; e.g, in the case of the operator J of total angular momentum 
the scalar symbol related to the matrix valued one (|3.6|) is given by 



j(p,x,n) = x x p + h + 1) n . (3.11) 

When applied to the principal symbol B (p,x), the relation ( p.lOj ) now associates 
with every quantum observable B of the type under consideration a classical observable 
bo(p, x, n), which is a function on the combined phase space ~R d x W 1 x S 2 . In particular, 
the Egorov theorem can be brought into the following form: According to ( |3.8| ) the matrix 
valued principal symbol of B(t) is of the same type as the right-hand side of (|2.16|) . Hence, 
due to the covariance property of the spin- Weyl correspondence with respect to SU(2), the 
associated scalar principal symbol bo{t)(p,x,n) = tr(Bo(t)(p,x)Aj(n)) reads 

b Q (t)(p,x,n) = bo(&(p,x),<p{d(p,x,t))n) = b (p(t) , x{t) , n(t)) . (3.12) 

The classical observable corresponding to the quantum observable B(t) at time t therefore 
follows from the respective classical observable at time zero by a purely classical time 
evolution on the combined phase space of translational and spin degrees of freedom. The 
combined flow on R d x M. d x S 2 , denoted as Y c *, is therefore given by the skew product, 

Y^p,x,n) := (&(p,x),<p{d(p,x,t))n) = (p(t) , x(t) , n(t)) , (3.13) 

and the scalar Egorov relation ( |3.12| ) reads 

b (t)(p,x,n) = b (Y*{p,x,n)) . (3.14) 



See, e.g., |2T[ for a general discussion of skew products. The spin part n i— > n{t) = 
ip(d(p, x, t))n of the combined classical dynamics takes place along the trajectories $'(p, x) 
of the translational motion and describes a precession of n(t) about the instantaneous axis 
G($'(p, x)), compare ( |2.18| ). This is a non-relativistic version of the Thomas precession 

of spin prrj] . 
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4 Quantum ergodicity 

Quantum ergodicity for Pauli Hamiltonians of the form ( p.2|) with spin 1/2 was derived in 
|7| under the condition that the underlying mixed classical/quantum system be ergodic. 
For this purpose, the two dynamics inherent in ( [3.8|) had to be combined into a single skew 
product flow Y* on R d x R d x SU(2) via 

y*((p, *),</) := (&(p,x),d(p,x,t)g) , (p,x,g) G R d x R d x SU(2) , (4.1) 

see also |10| for further details. Introducing Liouville measure (the microcanonical distri- 
bution) 



on the energy shell 



dfi E (p,x) :=—}—5(H (p,x)-E) dpdx (4.2) 

VOl \ I R 



VL E := {(p,x) G R d x R d \ H {p,x) = E) , (4.3) 



one observes that Y t leaves the product measure \i := |ijX fi H invariant, where /x# denotes 
the normalised Haar measure on SU(2). Now ergodicity of on x SU(2) with respect 
to \x means that for every integrable function F G ^(Qe x SU(2), d/x) the relation 

lim ± [ F(Y\t y, h)) dt = [ F(p, x, g) d//(p, x, g) (4.4) 



n B xsu(2) 



holds for /i-almost all initial conditions (£, y, h) G Q E x SU(2). Provided that ergodicity of 
Y** and certain properties of the principal symbol Hq of the quantum Hamiltonian hold, it 
was shown in J7| that in every sequence {ipk \ Ef, G I(E, h)} of orthonormal eigenspinors of 
the Hamiltonian Hp with eigenvalues Ek in an interval I(E, H) := [E — fiw, E + hu], which 
shall not contain any critical values of H , there exists a sub-sequence {ipk v \ Ek v G I(E, h)} 
of density one, i.e. 

#{„|£l6l(g,ft)} 
™ I £„ E /(£,»)} ' { ' 

such that for every quantum observable B with hermitian symbol B(p,x) and principal 
symbol B Q (p,x) the expectation values taken in the eigenstates {ipk v }f&N of Hp fulfill 

1 /" 1 
lim = o tr / B o(P> x ) dfj, E (p,x) =: -ti fi E (B ) . (4.6) 

^oo 2 2 

Moreover, the sub-sequence {V^jveN can be chosen independent of the observable _B. Here 
we remark that since the number Nj of eigenvalues G I(E, K) grows like 

(2 j + l)volfi B , s 
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(see fTO| , 0) in the semiclassical limit, both the numerator and the denominator in ( |4.5| ) 
become infinite as h — > 0. Thus, the limit v — > oo in (4.6) requires the semiclassical limit. 

We will now derive the analogous result for Pauli Hamiltonians with arbitrary spin 
under the condition that Y^, i.e. the purely classical skew product (|3.13|) , is ergodic on 
OgxS 2 with respect to the invariant measure /i c i := \ie x //g2. This means that for every 
integrable function / e L l (VL E x S 2 ,dyU d ) the relation 

1 ' T 



lim - / f(Y* 1 (£,y,m))dt= f(p, x, n) d/i cl (p, x, n) (4.8) 

i^oo 1 J Q Jn E xSV(2) 

holds for /i c i-almost all initial conditions (£, y, m) e Jig x S 2 . We remark that the ergodicity 
( [4. 3D of Y" f , which was required for quantum ergodicity in |7[, implies the relation ( |4.8|) , 
i.e. ergodicity of Y^, see appendix [FJ. Also notice that using property (iii) on page |I| the 
right-hand side of (|4.6|) can be rewritten as (J = =) 



-tr^(So) = / b (p,x,n) d// d (p, x, n) =: // c i(&o) . (4.9) 
1 Jn E xs 2 

In order to prove quantum ergodicity for arbitrary spin j along the lines of |22| , ||, [7|, 
we now analyse the quantity 

S 2 {E,h):=±- \{^k,B^ k )-^ cl (b )\ 2 , (4.10) 

1 E k &I{E,h) 

which is the variance of the expectation values of the quantum observable B about the 
classical mean of its scalar principal symbol bo. Introducing the bounded and self-adjoint 
auxiliary operator 

Br := i / ei^Be-^' dt - fi cl {b )± 2j+1 , (4.11) 
1 Jo 

the variance ( |4.1U| ) also reads 

S 2 (E,h) = ±- |(^,B T ^}| 2 <^- i^B 2 T ^ k ), (4.12) 

1 E k ai{E,h) 1 E k ai{E,h) 

where the estimate on the right-hand side follows from the Cauchy-Schwarz inequality. 
Quantum ergodicity is now implied, if this upper bound vanishes in the semiclassical limit. 
In order to show that this indeed happens, we employ the Szego limit formula 

fcF H (^,BV fc ) = ^— tr/i £ (B )=/i cl (6o) (4.13) 

for an observable B, which has been established in [[/]]. This formula is valid under certain 
technical assumptions on the symbol of the Hamiltonian Hp (see |7| for details), and 
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if the periodic orbits of the translational flow $* with non-trivial periods form a set of 
Liouville measure zero in Qe- Since shall be ergodic, which includes the ergodicity of 
the translational motion on the energy shell Qe, this condition is clearly fulfilled and we 
can apply (|4.13|) to (|4.12| ), yielding 



VtmS 2 (E,h) < — — -tr/i £ ((S Tj0 ) 2 ) = // d ((&r,o) 2 ) , (4.14) 



where on the rigth-hand side use has been made of the relation ( 2.13 ). Here Bta(p,x) 
and brflip, x, n) denote the matrix valued and scalar principal symbol, respectively, of the 
auxiliary operator ( 4.11]) 



In order to conclude the proof consider the scalar principal symbol 6r,o of Bt, which 
according to Q4.ll]) and (|3.14j) reads 



1 f 

b T ,o(p, x,n) = - J b {Y^(p, x, n)) dt - /i cl (b ) . (4.15) 

An application of the ergodicity ( |4.8| ) of Y* { to the principal symbol bo(p,x,n) of the 
observable B now reveals that &t,o(p, x , n ) vanishes as T — > oo for ^d-almost all (p, x, n) e 
Qe x S 2 . Since therefore the right-hand side of ( |4.14|) vanishes, this finally implies that 



lim S 2 {E, h) = . (4.16) 

According to a general construction in proofs of quantum ergodicity 0, the vanishing 
of S2(E,h) in the semiclassical limit implies the existence of a density-one sub-sequence 
{il>ku}veN C {ipk}ken along which one has 

\im(^ K ,B^) =/i c i(6 ) ■ (4.17) 

v—*oo 

Finally, a diagonal construction (see again 0, 0) yields a sub-sequence, still of density one, 
that is independent of the observable B such that ( f4.17p holds along this sub-sequence for 
all observables of the type considered here. 

Our conclusion therefore is that quantum ergodicity for Pauli Hamiltonians with arbi- 
trary spin j holds under the condition that the classical skew product Y^ of translational 
and spin dynamics is ergodic on flgxS 2 . We remark that ergodicity of Y£ implies ergodic- 
ity of the translational motion $' on the energy shell Qe- However, ergodicity of $* alone 
is not a sufficient condition for quantum ergodicity; see for a counter-example. 

For the purpose of an interpretation of the quantum ergodicity relation ( |4.17| ) for Pauli 
Hamiltonians with arbitrary spin, we now discuss its effect in terms of Wigner transforms 
of spinors ip G S'(R d ) ® C 2j+1 . Matrix valued Wi gner transforms are defined as 

W[ip](p,x) := [ e-^i){x-\y) ®^(x+\y) dy , (4.18) 

JR d 
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so that quantum expectation values of an observable B with symbol B(p, x) in states given 
by ip G L 2 (R d ) <g> C 2j ' +1 read 

(^Bi>) = —!—-[ [ ti (W[i/>](p,x)B(p,x)) dpdx. (4.19) 

After employing the relation ( |2.13| ) we obtain 

(ip, B^j) = 2 l + J~ I I I w[ip](p,x,n)b(p,x,n) dpdxdfi S 2(n) , (4.20) 
{2irh) cl J Rd J R d J S 2 

where 

w[ip)(p,x,ri) := tr (w[ip)(p,x) Aj(n)J (4.21) 

is introduced as a scalar Wigner transform of ^ £ 5'(M d ) (g>C 2j ' +1 . In this context, quantum 
ergodicity ( |4.17|) implies that along a density-one sub-sequence the scalar Wigner trans- 
forms {w['0Jfc v ]}i/eN °f eigenspinors of the quantum Hamiltonian become equidistributed on 
both the (translational) energy shell Qe and on the classical phase space S 2 of spin, i.e. 

lim j±±±w[M (p, x, n) = —L_ S ( Ho (p, x) - E) , (4.22) 



where the convergence has to be understood in the sense of distributions, i.e. after inte- 
gration with a suitable function on M. d x R d x S 2 as in ( 4.20Q . 
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A Explicit form of the kernel Aj(n) 

In this appendix we comment on some properties of the particular representation for the 
kernel Aj(n) that was chosen in |16j and give explicit examples for j = \ and j = 1. 

Expressing n 6 S 2 C I 3 in spherical coordinates, n = (sin 9 cos ip, sin 9 sin ip, cos 9) T , we 
denote the usual spherical harmonics by Yj m (n). Then the kernel Aj(n) can be expanded 
in the basis {Yj m (n) 1 1 G No, \m\ < /}, 

1=0 m=-l " J 



where for given j the C\ m are (2j + 1) x (2j + 1) matrices, see []T6 1 for details. Notice that 
this expansion terminates at I = 2j. 
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Thus, if one expands a general classical observable, given by some function b(n) G 
L 2 (S 2 ), in spherical harmonics and uses the spin-Weyl correspondence (|2.11|) in order to 



associate to it a spin-j quantum observable b, then its Weyl symbol (|2.12 ) will in general 



not coincide with the original b(n), but instead its expansion in spherical harmonics will 
be truncated at I = 2j. However, since the fundamental theory is quantum mechanics, one 
should adopt the reverse point of view: Starting from a spin-j quantum observable b one 
can assign to it its Weyl symbol b(n), given by (|2.12| ), as a classical observable. Then by 
( j2.11| ) this correspondence is one-to-one. 

Therefore, in a classical description of a spin j a general observable b(n) e L 2 (S 2 ) has 
an expansion in spherical harmonics Yj m (n) which terminates at I = 2j. Conversely, if the 
aim were a complete quantum description of a top whose classical dynamics is given by 
the Euler equations (instead of a classical description of, say, an electron with fixed spin 
|), one would have to perform the semiclassical limit h — > simultaneously with the limit 
of large spin, j — > oo, such that jh = 0(1), cf., e.g., |2"5| . 



As an explicit example we now briefly discuss some properties of Aj(n) in the repre- 
sentation chosen in |]TB| . For spin j = | we obtain 

A , (n)- r-(Y m (n)+Y 10 (n) -y/2 Y u (n) \_ x /7 
Al/2(nW7F l -V2Y u (n) Y 00 (n)-Y w (n))--^ + y* n ' (A ' 2) 

Since every observable b of a spin 1/2 is a hermitian 2x2 matrix, which can be represented 
as b = bo I2 + b ■ cr with bo and b e R 3 , the associated classical observable reads 



b(n) = tr(6 A 1/2 (n)) = b + V3b ■ n = b + 2y/ j(j + 1) b ■ n . (A.3) 



In particular, the spin operator S = ^cr is mapped to the symbol s = hy^3/An. This 
last statement generalises to arbitrary spin, i.e. S = |d7Tj(cr) is always mapped to s = 
fcyjij + 1) n - E.g., for j = 1 the inverse transformation (|2.11| ) reads 



S = 3 s Ai(n) d/xs2(n) 




y/2 \ / -iy/2 \ /2 0\\ T (A.4) 
y/2 V2 



y/2 

where now the kernel Ai(n) can be expressed as 




A x (n) = 1 3 + ^ (dTri(ir) • n) + ^f- {^{a) • n) 2 . (A.5) 



B Ergodicity of Y f and 

In this appendix we discuss the relation between the two skew products Y l ( |4.1|) and 
and their respective ergodic properties. 
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To this end we need a mapping of quantum spin precession on SU(2) to classical spin 
precession on S 2 . This is conveniently given by the Hopf map, one of whose possible 
realisations reads 

7i H : SU(2) — ► S 2 

tt (B.l) 
g i — >u { gt(Tgu, 

where u G C 2 with ||w||c 2 — 1 is arbitrary but fixed. We now want to relate the spin part 
g i — > d(p,x,t)g of the skew product flow Y l to the spin part n i— > <p(d(p, x, t))n of Y^. 
The definition of the universal covering map <p : SU(2) — > SO (3) implies that 

7T H (d(p, x, t)g) = p{d(p, x, t)) 7r H (g) , (B.2) 

so that the projection of the (quantum mechanical) time evolution in SU(2) to S 2 fulfills 
the classical equation of spin precession ( |2.18| ). This is a manifestation of the Ehrenfest 
theorem, since the Hopf map (p.2j ) yields the expectation value of (normalised) quantum 
spin in the state gu G C 2 , see ||. 

Therefore, defining ft H := id K d xR d®7r# : R d x W l x SU(2) -> R d x R d x S 2 , i.e. 
ir H (p,x,g) = (p,x,ir H (g)), we obtain the following commuting diagram: 

yt 

(p,x,g) ► (&(p,x),d(p,x,t)g) 



7TH 



yt 

(p,x,n) > (&(p,x),<p(d(p,x,t))n) 

Furthermore, to every (integrable) function / G L 1 {VLe x S 2 ,d/i c i) one can associate a 
function F G L l (VL E x SU(2), dfj) by defining F :— fo n H , i.e. F(p, x, g) = f(p, x, n H {g))- 
Provided now that Y t is ergodic on Qe x SU(2) with respect to fi = He x ^h, one observes 
that 

lim i f f(Y^y^ H {h))) dt = lim ± f F{Y\^y,h)) dt 

= / F(p,x,g) dfj,{p,x,g) (B.3) 

7n B xsu(2) 



/ f(p,x,n) dfi cl (p,x,n) 

Jn E xs 2 



for /i S -alomst all (£, y) G f2 s and /i#-almost all /i G SU(2). Since this implies the validity 
of (|B.3[) for /i S 2-almost m = 7r#(/i) G S 2 , the ergodicity of y* implies ergodicity of on 



fi^xS with respect to pL c \ = [Le X /xga. 
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